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ABSTRACT 

We present high spatial resolution ultraviolet and optical spectroscopy, ob- 
tained using the Space Telescope Imaging Spectrograph on board the Hubble 
Space Telescope, of nuclear structures within four 'warm' Ultraluminous Infrared 
Galaxies (ULIRGs). We find an AGN in at least three, and probably all four 
of our sample, hosted in a compact, optically luminous 'knot'. In three cases 
these knots were previously identified as a putative AGN nucleus from multi- 
band optical imaging. Three of the sample also harbor a starburst in one or 
more knots, suggesting that the optically luminous knots seen in local ULIRGs 
are the most likely sites of the dust-shrouded starburst and AGN activity that 
power the infrared emission. The four AGN have a diverse range of properties; 
two are classical narrow line AGN, one shows both broad and narrow lines and 
evidence for lines of sight from the narrow through to the broad line regions, 
and one is plausibly a FeLoBAL AGN. The probable presence in one object of 
an FeLoBAL AGN, which are extremely rare in the QSO population, supports 
the idea that LoBAL AGN may be youthful systems shrouded in gas and dust 
rather than AGN viewed along a certain line of sight. The three starbursts for 
which detailed constraints are possible show a smaller range in properties; all 
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three bursts are young with two having ages of ~ 4Myr and the third having 
an age of 20Myr, suggesting that ULIRGs undergo several bursts of star forma- 
tion during their hfetimes. None of the starbursts show evidence for Initial Mass 
Function slopes steeper than about 3.3. The metallicitics of the knots for which 
metallicities can be derived are all at least 1.5Zq. The properties of one further 
starburst knot are consistent with it being the forming core of an elliptical galaxy. 
Our results suggest that detailed studies of the 'knots' seen in ULIRGs can give 
important insights into the most violent starburst and AGN activity at both low 
and high redshift. 

Subject headings: galaxies: active — galaxies: Seyfert — galaxies: starburst — 
galaxies: interactions — galaxies: nuclei — ultraviolet: galaxies 

1. Introduction 

Ultraluminous Infrared Galaxies (ULIRGs, those objects with l-1000//m luminosities 

in excess of lO^^L©) were first discovered by Rieke & Low (1972), and have been intensely 
studied since the discovery by the Infrared Astronomical Satellite (IRAS) of significant pop- 
ulations of ULIRGs in the local Universe (Soifcr ct al 1984). Unveiling the power source 
for the infrared (IR) emission in ULIRGs, thought to be some combination of starburst and 
AGN activity triggered by ongoing interactions (see Sanders & Mirabel (1996) for a review), 
is important for two reasons. Firstly, local ULIRGs, although rare, represent the closest 
examples of the most extreme cases of star formation and black hole accretion in galaxies, 
and can be used to study in detail how such processes drive galaxy evolution. Secondly, the 
ULIRG luminosity function shows strong evolution with redshift (Saunders et al 1990) and 
deep sub-mm surveys (Hughes et al 1998; Eales et al 2000; Scott et al 2002; Borys et al 2003) 
have found that systems with ULIRG-like luminosities are much more numerous ai z > 1 
than locally. A detailed understanding of local ULIRGs can therefore give insights into the 
power source (Rowan-Robinson 2000; Verma et al 2002; Farrah et al 2002b, 2004) and the 
role of interactions (Farrah et al 2002a; Chapman et al 2003; Conselice et al 2003a,b) in 
distant IR-luminous sources, and the impact of these sources on the cosmic history of star 
formation. 

After just over three decades of study, a consensus has begun to emerge on the origin 
of the IR emission in ULIRGs. Early studies using optical spectroscopy showed that some 
ULIRG spectra resembled those of starburst galaxies (Joseph & Wright 1985), whereas others 
contained emission lines characteristic of Seyferts (Sanders et al 1988a). Radio observations 
showed direct evidence for dust shrouded starbursts in some ULIRGs (Condon et al 1991; 
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Smith et al 1998) and AGN cores in others (Lonsdale et al 2003; Nagar et al 2003). Papers 
based on mid-infrared spectroscopy suggested that most ULIRGs (~ 80%) are powered by 
starbursts (Rigopoulou et al 1999; Genzel et al 1998), but that at least half of local ULIRGs 
contain both a starburst and an AGN. Results from modelling the 1-1000/im spectral energy 
distributions of local ULIRGs (Farrah et al 2003) show that all ULIRGs contain a dusty 
starburst, but that only about half contain an AGN that contributes significantly to the 
total IR emission. These results are consistent with X-ray observations of small samples of 
ULIRGs, which indicate that the bolometric emission from most ULIRGs is dominated by 
starburst activity rather than an AGN (Pranceschini et al 2003; Ptak et al 2003). There is 
also evidence that the presence of an AGN in ULIRGs depends on the total IR luminosity; 
optical and near-IR spectroscopic surveys (Veilleux et al 1995; Veilleux, Sanders & Kim 
1999) show a sharp rise in the fraction of ULIRGs with AGN signatures with increasing IR 
luminosity, and mid-IR spectroscopy from ISO (Tran et al 2001) showed that ULIRGs with 
Lir < 1O^^-''L0 were mostly starburst dominated, whereas ULIRGs with iLj^ > 10 ' Lq were 
more likely to contain an AGN. 

Although we now have a much clearer picture of the role of starbursts and AGN in 
ULIRGs, few studies have, to date, correlated starburst and AGN activity in ULIRGs with 
morphological features (as opposed to correlating the power source in ULIRGs with broad 
scale morphologies, see e.g. Veilleux, Kim, & Sanders (2002)), due mainly to the poor spa- 
tial resolution of current mid- and far-IR instruments when compared to the sub-arcsecond 
scale structures seen in UV/optical images of local ULIRGs. Correlating power source with 
morphological features in ULIRGs however has the potential to be a powerful probe of their 
nature and evolution; most ULIRGs display a plethora of tails and loops, and compact 
'knots' in the nuclear regions (Surace et al 1998; Farrah et al 2001; Bushouse et al 2002), 
and the 'knots' in particular are suggested by some observations as possible sites for the 
dust-shrouded starburst and AGN activity in ULIRGs (Surace et al 1998; Soifer et al 2000). 

A subset of the ULIRG population that is well suited to exploring links between mor- 
phology and power source are those ULIRGs with 'warm' infrared colors (/25//6O > 0.2)^. 
This color selection favors ULIRGs that contain dust at temperatures > 50K, which is most 
plausibly attributed to an AGN rather than a starburst. Coupled with the apparent ubiquity 
of starbursts in ULIRGs (Farrah et al 2003), this color selection is hkely to include those 
ULIRGs that contain both a starburst and an AGN, where the interplay between these two 
processes can be studied in detail. In this paper, we use the Space Telescope Imaging Spec- 
trograph (STIS) on board the Hubble Space Telescope (HST) to obtain long-slit ultraviolet 



^The quantities /12, /25, /eo & /loo represent the non color corrected IRAS flux densities in units of 
Janskys at 12/ini 25/um 60/xm & lOO/xm respectively 
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and optical spectra of several 'knots' within the nuclear regions of four 'warm' ULIRGs. 
Observations and analysis are described in §2 and results are given in §3. Our conclusions 
are presented in §4. We assume f2 = 1, A = 0.7 and Hq = 65 km s~^ Mpc~^. Luminosities 
are quoted in units of bolometric solar luminosities, where Lq — 3.826 x 10^^ Watts. 

2. Observations &; analysis 

We selected for observation ULIRGS from the sample with 'warm' mid-infrared colors 
imaged in the optical using the Wide Field Planetary Camera 2 (WFPC2) on board HST 
by Surace et al (1998) and with adaptive optics by Surace & Sanders (1999). This imaging 
demonstrated the presence of compact (<100 pc) high surface brightness features in the 
galaxy cores. Furthermore, these 'warm' systems show evidence for Scyfcrt activity, which 
led Sanders et al (1988a) to postulate that these were objects transitioning between the 
general ULIRG population and optically selected quasars. Increasing central concentration 
in the infrared, as well as optical colors and luminosities, led Surace et al (1998) and Surace 
& Sanders (1999) to further identify some of these knots as possible active nuclei. 

The four objects that were observed (Table 1) span a range in properties, including 
both red (F12071-0444) and blue (F01003-2238) colors, simple (F01003-2238) and com- 
plex (F15206+3342) nuclear morphologies, single (F01003-2238, F12071-0444) and double 
(F05189-2524) nucleus systems, and HII, Sy 2, and Sy 1.5 spectral types. The data were 
taken in cycle 8 (proposal ID 8190) with STIS in longslit spectroscopy mode on board the 
HST. Three spectra were taken of each object, with the 52" x 0.2" slit and the G140L, 
G470L and G750L gratings, providing coverage over the wavelength ranges lOOOA - 1700A 
and 3000A - lOOOOA in the observed frame. Spatial resolutions are 0.024 arcseconds pixel"^ 
for the UV spectra, and 0.05 arcseconds pixel"^ for the optical spectra. Spectral resolutions 
are a FWHM of 1.5 pixels at 1500A for the G140L grating, 1.3 pixels at 3200A for the G430L 
grating, and 1.6 pixels at 7000A for the G750L grating. 

In all cases the slit was centered on the putative nucleus whose luminosity and color were 
consistent with an obscured AGN (Surace et al 1998), and oriented so it also covered one 
or more off-nuclear 'knots', if such knots were present. In Figure 1 we present the WFPC2 
images of our sample, with the STIS slit positions superimposed. In Table 2 we list the 
magnitudes and radii of these knots, together with the stellar masses inferred by Surace et al 
(1998) under the assumption that all of the B and I band flux is stellar in origin. The total 
exposure times were 5500s for the G140L spectra, 2700s for the G470L spectra, and 1750s 
for the G750L spectra. Each observation was split into three or more exposures dithered 
with integer pixel shifts to allow the removal of cosmic rays and bad pixels. 
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Individual exposures were reduced using the automated calstis pipeline at the STScI and 
the latest calibration files available as of January 2004. The exposures were then combined 
into a single image using the Image Reduction and Analysis Facility (IRAF) task ocrreject. 
Spectra were extracted using apertures corresponding to the same physical size in the UV and 
optical for each object. For F01003-2238, F12071-0444 and F15206+3342, apertures of width 
3 pixels for the optical spectra and 7 pixels for the UV spectra were used, corresponding to a 
size of approximately 350pc. For F05189-2524, where the two knots covered by the STIS slit 
lie ~ SOOpc from each other, apertures of width 5 pixels for the optical spectra and 11 pixels 
for the UV spectrum were used, corresponding to a smaller physical size of approximately 
220pc. 

Line profiles and errors were measured using the spec fit package within IRAF. This 
package can fit a wide variety of emission profiles, absorption profiles and continuum shapes 
to an input spectrum, and can disentangle multiple blended lines. The fitting itself was 
carried out via minimization, using a combination of Levenberg-Markwardt and Simplex 
algorithms. Input parameters for specfit were determined using the IRAF task splot. 

To deredden the UV spectra the color excess of the ionized gas, E{B — V)gas was derived 
from the observed and (assumed) intrinsic Balmer decrements: 



£(B-VU^ 0.935 xtajl^ (1) 

where intrinsic Balmer decrements of 2.85 and 3.1 were assumed, corresponding to heating 
by stars and an AGN respectively. This color excess was then converted to the reddening of 
the stellar populations via (Calzetti 2001): 



EiB - VUars = OAAE{B - V)gas (2) 

The extinction at some wavelength A, A^, can then be computed using the prescription given 
by Leitherer et al (2002): 



Ax = E{B - V)stars X 

. .rr^ 0.671 9.218 X 10-3 2.62 x lO'^N 



This formula is valid over the wavelength range 0.097 < X{fj,m) < 0.18. The optical spectra 



-6- 



were dereddened using the prescription given by Cardelli et al (1989)^. 

Line Full Widths at Half Maximum (FWHM) were corrected for the effects of finite 
instrument resolution using the quadrature method, where: 

FWRM-"^^^, = FWHMl^ - FWHM^^ (4) 

This method assumes that the real and instrument profiles are Gaussians, and therefore gives 
line widths that are systematically too high for profiles that are more peaky than Gaussians. 
As such peaky line profiles are known to be present in AGN (Veilleux 1991), the discussion 
of line widths in this paper should be regarded with caution. 

3. Results 

The UV and optical spectra for each ULIRG are presented in Figures 2, 3, 4, 5, 6, 7 
& 8. The UV and optical emission line fluxes are presented in Tables 3 & 4 respectively. 
Widths of selected emission lines are given in Table 5. 

All of the spectra show at least two emission lines, while a few also show absorption 
lines. These lines are a powerful tool in determining physical conditions in the ionized gas 

in our sample. In the optical, diagnostics using the stronger emission lines can be used to 
constrain the temperature, density and kinematics of the hot gas in galaxies. In this paper 
we have used the three line ratio diagnostics presented by Veilleux & Ostcrbrock (1987) 
as the primary method to distinguish between HII and AGN spectra; [OIII] A5007/H/9 vs. 
[Nil] A6583/Ha, [OIII] A5007/H/? vs. [SII] AA6716, 6731/Ha, and [OIII] A5007/H/3 vs. [01] 
A6300/Hq;. For those spectra classified as AGN we adopt the definition that a LINER is 
an AGN with log([OIII] A5007/H/3) < 0.477. For those spectra classified as HII regions we 
have used the theoretical upper limits on these line ratio diagnostics for both continuous star 
formation (Kewley et al 2001a) and instantaneous bursts (Dopita et al 2000). All of these 
diagnostics are presented in Table 6, and plotted in Figure 9. We supplemented these with 
two additional diagnostics; [OIII] A5007/H/? vs. [Oil] A3727/[OIII] A5007 (Baldwin, Phillips 
& Terlevich 1981) and [OIII] A5007/H/3 vs. [Oil] AA7320, 7330/Ha (Osterbrock, Tran & 
Veilleux 1992). 



^This paper also presents a prescription for dereddening UV spectra at wavelengths as small as O.l/Ltm, 
however their far-UV segment is derived using limited data, and therefore we do not adopt their prescription 
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Spectroscopy in the ultraviolet, despite suffering significant extinction uncertainties, is 

the best way of directly studying the hottest, youngest stars, as at wavelengths longward 
of 3000A features from these stars are weak, and are often blended with nebular emission 
lines. The shapes of the stronger UV lines can be used to distinguish between starburst 
and AGN activity and shocks vs. photoionization, and determine parameters for the stellar 
initial mass function (IMF). Here we have analysed the UV spectra of our sample using both 
UV line ratio diagnostics, and modelling using the spectral synthesis package Starburst99 
v4.0 (Leitherer et al 1999). 

In the optical, star formation rates are commonly estimated using the fluxes of recom- 
bination lines, particularly Ha, calibrated using spectral synthesis models (see Kennicutt 
(1998) for a review). Though the scatter in the cahbration can be up to 30%, due to differ- 
ences between spectral synthesis models, this method has the advantages of high sensitivity, 
and direct coupling between the Ha flux and the (massive) star formation rate. The pri- 
mary disadvantage, for the purposes of our study, is that the Ha flux from stars must be 
disentangled from the nebular [NII]AA6548, 6583 lines, and from any broad Ha line arising 
from an AGN. In most cases it was not possible to reliably isolate the Ha flux arising from 
star formation (with errors of a factor of two or more, see §§3.1 - 3.4 for more details); we 
have therefore, for consistency, estimated optical star formation rates for all objects using 
the extinction corrected luminosity of the [OH] A3727 line, and the prescription presented 
by Kennicutt (1998): 

SFR = 1.4 ± 0.4 X 10-^^L([O//]) (5) 

where the SFR is in units of Mq yr~^ and L[OII] is in units of ergs s""^. The luminosity 
of [OH] A3727 is sensitive to the metallicity and ionization state of the gas in the narrow 
line region, and the calibration for [OH] A3727 as a star formation rate indicator is derived 
empirically through relationships derived for Ha. The errors on the star formation rates 
from Equation 5 are therefore at least 30%, and these star formation rates should therefore 
be regarded with caution. We also emphasize that, for those 'knots' that contain an AGN, 
the star formation rates derived using Equation 5 should be regarded as upper limits. 

Where possible, we estimated electron densities and temperatures using the IRAF task 
temden. The commonly used line ratios for estimating electron temperatures are [OIII] 
AA4959, 5007 / [OIII] A4363 and [NII]A5755/[NII]A6583, however [OIII] A4363 and [NII]A5755 
are not detected in any galaxy of our sample. Electron temperatures were therefore estimated 
using the line ratio [OH] A3727 / [OH] A7325. To estimate electron densities we used the 
[SII] A6716 / [SII] A6731 line ratio, though this was only possible in a few cases due to the 
blending of the two [SII] lines. 
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3.1. IRAS F01003-2238 

3.1.1. Background 

This object contains a single nucleus embedded in a small halo of extended emission, with 
a series of small blue knothke features extending 1" to the north-west (Surace et al 1998). 
A serendipitous observation of this object by ROSAT did not detect any X-ray emission, 
and gave a la upper limit to the soft X-ray luminosity of 2.14 x 10^° ergs s~^ (Stevens & 
Strickland 1998). Previous efforts to identify the power source in this galaxy have produced 
conflicting results. An early study by Armus, Heckman & Miley (1988), using ground-based 
spectroscopy of the whole galaxy, identified several features characteristic of Wolf-Rayet 
stars (Hell A4686, Hel A5876, and tentatively NIII AAA4634, 4640, 4642 GUI A5696, NIV 
A5737, and several features around 7900A, though see also Conti (1991)), and concluded 
that this galaxy contains approximately 10^ WN-type Wolf-Rayet stars, consistent with 
results from 2-D spectroscopy (Wilman, Crawford & Abraham 1999). Surace et al (1998) 
however postulate that the central nucleus is an AGN, based on its position in optical colour- 
colour plots, and mid-infrared spectroscopy suggests that the IR emission is not dominated 
by a starburst (Tran et al 2001). Near-infrared spectroscopy (Veilleux, Sanders & Kim 1997) 
however shows no evidence for a hidden broad line region. 

3.1.2. Emission & absorption line properties 

The data show a single, bright spectrum in the UV and optical bands. In the UV 
spectrum Lya emission is present, together with weaker emission lines of NV A1240, SilV 
AA1394, 1403, and NIV] A1488. There are also a number of absorption systems, and weak, 
but distinct PGygni profiles on the SilV doublet and NV lines, which imply winds from O 
supergiants. The SilV doublet is commonly seen in O star photospheres, however the NIV] 
A1488 fine is only seen in winds from W-R stars. We also clearly detect the Wolf-Rayet 
star signature triplet NIIIAAA4634, 4640, 4642 (which may be blended with GIV A4648). We 
do not however detect the 'bump' between 7875A and 7965A claimed by Armus, Heckman 
& Miley (1988). This could indicate either that this feature does not exist, or that it is 
produced outside the nuclear region. The B band spectrum contains the [NeV]A3426 line, 
which is thought to originate in AGN. The emission lines in the UV and optical spectra are 
spatially symmetric with respect to the continuum, with extents of 550-600pc for the Ly« 
and stronger optical lines, and ~ 280pc for the high ionization UV lines (for F01003-2238 
one pixel in the UV frame corresponds to ~55pc and 1 pixel in the optical corresponds to 
~110pc). The fine profiles of H/5 and [GUI] are also symmetric, and show no blue wings. 



-9- 



This is at odds with ground-based spectra (Armus, Heckman & Miley 1988; Veilleux, Sanders 
& Kim 1999; Lipari et al 2003) which show a double peaked profile on H/3 and [OIII], with 
an inferred outflow velocity of ~ 1000km s^^. As our spectra are of comparable spectral 
resolution, we would expect to sec these wings, if they are present. The implications from 
this are discussed in §3.1.4. The ground-based spectra also show a fiat continuum, whereas 
our optical spectra are very blue; this difference can be ascribed to our spectra targeting 
only the active nuclear regions rather than the whole galaxy. 

Three lines in the 1 band spectrum, at (observed-frame) 6416A, 6535A and 694lA, 
proved difficult to identify reliably. The line at 694lA has no plausible IDs at the sys- 
temic redshift (z=0.117). This line could be [FeVll]A6087, but this implies a velocity shift 
of 7000 ± 1500km relative to systemic (or a redshift of z=0.141). This shift is, to 
our knowledge, at least 3000km s"^ larger than the shifts seen for this hne in any other 
galaxy. Adopting this same shift, then the features at 6416A and 65 35 A would be identi- 
fied as [FcVI]A5632 and [FcVII]A5721 respectively. In this circumstance the required ratio 
[Fe\///]A6087/[FeV//]A5721 > 1.0 is also satisfied. Equally plausible IDs for the features 
at 6416A and 6535A however are CIIIA5696 and CIVA5808 respectively, both signatures of 
WC type Wolf-Rayet stars. In this case the implied velocity shift is a perfectly reasonable 
~ 1300km s~^ (or a redshift of z=0.125). Finally, plausible IDs for the features at 6416A 
and 6535A at the systemic redshift are [NII]A5756 and HeIA5876. 

These three lines at 6416A, 6535A and 694lA present a conundrum. All three can be 
identified as high excitation iron lines at the same redshift, but the implied velocity shift 
is barely believable. Two of the three features can also be identified as lines characteristic 
of WC type W-R stars, which is consistent with the other W-R star signatures seen in the 
UV and optical spectra. These same two lines can also be identified as commonly seen 
emission lines at the systemic redshift. Given this uncertainty, any conclusions we draw are 
uncertain. Barring a plausible alternative for the feature at 694lA we tentatively conclude 
that this feature is [FeVII]A6087. This implies that the other two feature may be due at 
least in part to high excitation iron lines, but the presence of NIIIAAA4634, 4640, 4642 also 
implies that they could be due to W-R stars. Without higher spectral resolution data, we 
can only adopt the most conservative conclusions; the line at 6535A we conclude is a blend 
of [FeVII]A5721, CIVA5808 and HeIA5876, and the feature at 6416A we conclude is a blend 
of [FeVI]A5632, CIIIA5696 and [NII]A5756. 

A careful search for the more common stellar absorption features (Ca(H)A3934 and 
Ca(K)A3968; the 'G' band feature at 4250A due to CN and CH; Mg I b A5176,5200; Na 
IDA5896, dTiO AAA6180,7100,7700; Call A8542) revealed only weak G band absorption, 
and a possible absorption feature from Call A8542. Their weakness or absence in these 
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spectra is not unexpected; our observations only include the small part of the host galaxy 
that lies along the line of sight to the nucleus, and these features are known to be diluted 
by either an AGN power-law continuum, or the continuum from hot young stars. 

3.1.3. Wolf-Rayet stars 

The presence of W-R star spectral features allows us to set hmits on the stellar popu- 
lations within the nuclear regions (see e.g. Schaerer & Vacca (1998)). Armus, Heckman & 
Miley (1988) postulate, based on the comparable strength of the NIII AAA4634, 4640, 4642 
lines and the Hell A4686 line, and the absence of NV A4604 and CIV AAA4648, 5801, 5812, 
that the W-R stars in F01003-2238 arc mainly WN stars, with relatively few WC stars. 
Our spectra are consistent with the presence of a large population of WN stars, but for the 
reasons given in §3.1.2 we cannot draw reliable conclusions about the presence or otherwise 
of WC stars solely on the basis of hne identifications. Our spectra however do not show the 
oxygen hues that are characteristic of WO stars (OIV A3400, OVI AA3811, 3834, Barlow & 
Hummer (1982)). We therefore conclude that at least some of the WR stars in the nucleus 
of F01003-2238 are of the WN subtype, with some unquantified fraction being WC stars. 
Concerning ourselves solely with the WN stars, then identifying the subtype of WN star is 
difficult. Again based on the relative strengths of the Hell A4686 line and the NIII 'bump' 
at 4600A, we postulate that the WN stars in F01003-2238 arc probably late type (i.e. WNL 
stars), in agreement with Armus, Heckman & Miley (1988), though higher spectral resolu- 
tion data to resolve the NIII AAA4634, 4640, 4642 lines would be required to confirm or refute 
this. 

With an identified WR star subtype, the number of WR stars can be estimated, as can 
the ratio of the number of WR stars to the number of O stars (strictly speaking the number 
of O stars of sufficient mass to evolve into WR stars, typically > 3OM0). We here assume, 
solely for the purposes of deriving Nwr/Nq, that all the W-R stars are WNL stars. In 
the solar neighborhood, Nwr/Nq ~ 0.1, a value consistent with continuous star formation. 
Higher values of Nwr/Nq can only be accommodated if the star formation happens in short 
duration bursts, or if the metal abundance is above Solar. The number of WNL stars can be 
estimated using the flux of the Hell A4686 line (Vacca & Conti 1992), assuming that all of 
the Hell A4686 flux is due to WNL stars. Taking the Hell A4686 luminosity of a single WNL 
star to be 1.7 x lO^^ergs s^^, then for our spectra this gives 1.2 x 10^ or 1.4 x 10^ WNL stars 
for intrinsic Balmer decrements of 2.85 and 3.1 respectively. Computing the ratio of WNL 
stars to O stars is more complex. The prescription given by Vacca & Conti (1992) for local 
W-R galaxies assumes pure case B recombination, and a single effective temperature and 
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electron density. The maximum value of Lo6s(A4686)/Lo6s(H/3) from Vacca & Conti (1992) is 
0.21, corresponding to excitation purely by WNL stars, whereas our derived value is 0.24. We 
therefore conclude that our initial assumption is not valid, and that F01003-2238 contains 
a significant number of WC stars in addition to the WN stars. Given the uncertainties in 
reddening introduced because the Ha and [N1I]AA6548,6583 lines could not be resolved, and 
the probable presence of an AGN (§3.1.4), reliable estimates of the ratio of W-R stars to 
O stars are not possible without higher spectral resolution data. Nevertheless, our results 
imply at least 10^ W-R stars in the nucleus of F01003-2238, and a value of Nwr/Nq that 
is substantially higher than in the Solar neighborhood; this is supportive of a 'burst' rather 
than continuous star formation. 



3.1.4- Analysis 

The I-band spectrum lacks the spectral resolution to rehably deblend the Ha and [Nil] 
lines, with variations in the derived Ha flux of up to a factor of 3 depending on the ini- 
tial line parameters input to Specfit. We therefore adopted two canonical values for the 
[Nil] X6583/Ha hne ratio, [Nil] A6583///a = 1.0 (corresponding to an AGN) and [Nil] 
A6583/ifQ; = 0.25 (corresponding to an HII spectrum), and derived line ratio classifications. 
In both cases, the diagnostics, presented in Table 6 and Figure 9, are ambiguous. Other diag- 
nostics do not clarify matters; the [0111] X5007/H(3 vs. [Oil] A3727/[0111] A5007 diagnostic 
classifies the spectrum as an Hll region, and the [0111] X5007 /H/3 vs. [Oil] AA7320, 7330/ Ha 
diagnostic classifies the spectrum as an AGN, irrespective of the assumed intrinsic Balmer 
decrement. The extinction corrected star formation rate derived using Equation 5 is 6±2Mq 
yr~^ and 100 ± 30Mq yr~^ for assumed intrinsic Balmer decrements of 2.85 and 3.1 respec- 
tively. 

The [SII] A6716 and A6731 hnes could not be deblended in this object, hence we only 
present an estimate of the electron temperature in the narrow line region. For an assumed 
intrinsic Balmer decrement of 2.85 the [Oil] A3727 / [Oil] A7325 ratio is 8.11, whereas 
for an intrinsic Balmer decrement 3.1 the ratio is 1.7. For a ratio of 1.7, the electron 
temperatures are 18, 700K, 6, 600K, and 5, 300K for assumed electron densities of Ue = 
2 X lO^cm"^, He = 1 X lO^cm"^ and rie = 2 x lO^cm"^ respectively. For a ratio of 8.11, 
the electron temperatures are lower, at 12000K, 8550K and 6500K for assumed densities of 
He = 5 X lO^cm"^, He = 1 X lO^cm"^, and rig = 2 x lO^cm^'^ respectively. This indicates 
that the narrow line region in F01003-2238 does not contain particularly dense gas, and 
that shock heating does not play a significant role. The detection of the [Oil] A7325 line is 
however marginal, at ~ 3a; these results should therefore be regarded with caution. 
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The clear signatures of massive star formation in the UV spectrum allow us to examine 
the starburst activity in F01003-2238 in detail. We have modeled the UV spectrum using 
the Starburst99 v4.0 spectral synthesis code (Leitherer et al 1999) using model parameters 
spanning the following ranges; the mode of star formation (i.e. continuous or a 'burst'), 
a starburst age of between IMyr and lOOMyr, an initial mass function slope of between 2 
and 30, and a lower stellar mass limit of between O.IM© and 8Mq. The following analysis 
assumes that the star formation in F01003-2238 can be described by a Salpeter IMF, and 
that the stellar evolutionary tracks used by Starburst99 are apphcable. The dereddened UV 
spectrum, together with the best-fit Starburst99 model, are plotted in Figure 11. The W-R 
star lines, together with the strong PCygni profile on the NV fine and SilV doublet, provide 
strong constraints on the starburst age, limiting the age to 3-4Myr (corresponding to the 
time taken for a 40Mq O star to evolve into a W-R star). A 'burst' of star formation is also 
favored over continuous star formation. The youth of the starburst makes constraining the 
lower stellar mass limits difficult; we found that the data were consistent across the whole 
O.IM0 to 8M0 range. Constraints can however be set on the slope of the IMF: the PCygni 
profiles, together with the CIII A1175 absorption feature, require an IMF slope of 3.3 or less. 
IMF slopes steeper than about 3.3 predict weak or no PCygni profiles or absorption troughs 
in young starbursts, and can be excluded. There are however two significant differences 
between the observed spectrum and the best-fit model. Firstly, the Lya emission seen in the 
spectrum is not in the model, as Starburst99 models younger than about lOMyr in general 
show Lya in absorption. As the Lya line is broad, with a corrected FWHM of ~ 2100km s~^, 
we attribute this line to an AGN rather than the starburst (though we note that this width is 
only slightly over the 'official', but entirely arbitrary lower limit for line widths in broad line 
regions of 2000km s~^). A corollary to this is that, if high excitation iron fines are present in 
the I band spectrum, then this would also indicate the presence of an AGN (a full discussion 
of possible origins of the 'coronal' iron lines is presented in §3.2). Secondly, the SiIIA1260 
and CIIA1335 fines seen in the model are not seen in the spectrum. We attribute this to a 
combination of lack of instrumental resolution, and dilution from the AGN continuum. 

The non-detection of blue wings on the H/3 and [OIII] lines is intriguing, given that 
such wings are clearly detected in lower resolution ground based spectra. The simplest 
conclusion would be that these blue wings do not exist, however they were detected in 
two independent ground-based measurements (Armus, Heckman & Milcy 1988; Lipari et al 
2003). This apparent conflict is best explained in terms of different apertures; the STIS 
spectra cover only the central 350pc, whereas the ground-based spectra cover 2kpc or more. 
If the blue wings seen in the ground-based spectra are interpreted as evidence for an outffow, 
then two possibilities exist. The first is that the outfiow does not originate in the nucleus. 
We cannot discount this possibility, however we consider it unlikely. The outfiow found in the 
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ground-based spectra is almost certainly driven by a starburst or AGN (the other possibihty 
being that it is a dynamical relic from the merger); we expect starbursts and AGN to be 
triggered in the nuclei of ULIRGs, and FO 1003-2238 is a single nucleus system. The second 
possibility is that this wind was generated some time ago, and that the STIS aperture only 
covers a small part of the wind. Taking a mean velocity of 1000km s~^ and assuming a 
distance of 500pc then this implies that the wind was initiated around 500,000 years ago. 



3.2. IRAS F05189-2524 

3.2.1. Background 

This object is classified as a late-stage merger from WFPC2 and ground-based imaging 
(Surace et al 1998; Veilleux, Kim, & Sanders 2002). Ground based observations show a Sy2 
spectrum (Veilleux et al 1995), and what appears to be a single, very red nucleus bisected by a 
dust lane (Surace Sz Sanders 1999; Scoville et al 2000). The presence of an obscured AGN was 
confirmed both by the observations of broad lines in polarized flux (Young et al 1996), and in 
direct light (Veilleux, Sanders & Kim 1999). Later observations (Dudley 1999) observed the 
11.3;um dust feature, suggesting that this object also contains an obscured starburst. The 2- 
lOKcV X-ray spectrum (Risaliti et al 2000) is well fit by a two-component model, consisting 
of a power law with F = 1.89^q;34, absorbed by a column density of Nh = 4.7 x 10^^^^;^ 
cm~^, and a thermal component with kT = 0.881q 35 keV. The power law component can be 
interpreted as arising from a Compton thin AGN, with the thermal component either due 
to the AGN or to starburst activity. Results from fitting the 1 — 1000//m spectral energy 
distribution of this object (Farrah et al 2003) found that both a starburst and an AGN 
were required to explain the infrared emission, with the starburst as the major contributor. 
Conversely, recent X-ray observations (Imanishi et al 2003) suggest that the starburst is 
not energetically dominant, with the AGN as the probable source of most of the bolometric 
emission. 

As shown in Figure 1 the STIS slit was placed across both bright regions in the center of 
F05189-2524. Consequently, two spectra are visible in the STIS image, which are discussed 
in the following subsections. 
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3.2.2. Knot one 

Emission line properties 

This knot is the brighter of the two in both the UV and the optical. In the UV several strong 
emission lines can be seen, including Lya and CIVA1549. Several weaker emission lines are 
also present, including NIV] A1488. This hne is not seen in AGN, including composite 
QSO spectra (Vanden Berk et al 2001) and is thought to originate in winds from Wolf- 
Rayct stars. None of the UV lines show PCygni profiles or evidence for absorption. The 
Lja line is spatially extended, and asymmetric with respect to the continuum, whereas the 
spatially compact higher ionization lines, including NVA1240, show no asymmetry (Figure 
10). Asymmetry is also seen in the Ha, H/9 and [0111] lines. The spatial extent of the Lya, 
Ha, H/? and [OIII] hues is approximately 450 - 500pc, whereas for the higher ionization UV 
lines the spatial extent is 250 - 300pc. The asymmetric region on the Lya and optical hues 
is ~ 130pc across (for F05189-2524 the STIS pixel scale is ~ 22 pc pixel"^ in the UV and 
--^ 50 pc pixel"^ in the optical). We interpret this region as evidence for asymmetry in the 
ionized gas in the nuclear regions, possibly a sign of a late stage merger. In the UV the 
line widths are broad, with a corrected Lya FWHM of 2175km s^^. In the optical both a 
narrow and broad component can be seen on the HjS line. None of the optical or UV lines 
show blue wings, though the broad components of H/3 and [OIII]A5007 show a systematic 
blueshift compared to the other optical lines, corresponding to a velocity of 300-600 km s~^ 
for H/3 and 900-1200 km s"^ for [OIII]A5007. The widths of the 'narrow' hues are broad, 
with corrected FWHMs of order 1000 km s~^. No W-R star features are visible in the optical 
spectra. 

In addition to the standard emission lines (Figure 3 and Table 4), three features are 
present in the B band spectrum between 3050A and 3350A that are not commonly seen in 
starburst or AGN spectra. The features at ~ 3080A are most likely to be Fell multiplets; as 
they lie at the extreme edge of the spectrum we do not consider them further. The line at 
3250A we identify as the OIII A3133 resonance-fluorescence line. Our reasoning is as follows. 
Adopting a redshift based on the Lya line gives a rest-frame wavelength for this feature of 
313lA, which is the same within the errors as the wavelength of the strongest of the Bowen 
resonance-fluorescence lines of OIII, at 3133A . The OIII resonance-fluorescence lines are 
created in nebulae via excitation from Hell Lya photons with rest-frame wavelength 303. 8A 
(Osterbrock 1989), so there is a physical motivation for adopting a redshift based on the 
UV lines rather than the optical lines. The OIII A3133 line is also seen in composite QSO 
spectra (Vanden Berk et al 2001). Finally, we identify the line at 3333A as being the Hel 
3 189 A triplet. 
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High excitation iron lines 

A number of high excitation iron hnes, hsted in Table 7, are present in the optical spectra. 
These lines show a systematic blueshift, corresponding to a velocity of ~ 1200km s^^, and 
include the very high ionization states [FeX] A6375, [FeXI] A7892 and, remarkably, [FeXIV] 
A5303 (which may be blended with [CaV] A5309). The [FeX]- [FeXIV] lines are usually only 
seen in spectra of local supernova remnants, and in the Solar corona (hence their description 
as 'coronal' lines), and are rare in extragalactic objects. To our knowledge, only two other 
galaxy spectra show a direct detection of [FeXIV]; III ZW 77 (Osterbrock 1981) and MCG- 
6-30-15 (Reynolds et al 1997). 

Two possible origins have been proposed for these coronal iron lines: (1) they originate 
in a 'Coronal Line Region' (CLR) intermediate in distance between the broad and narrow 
line regions (Penston et al 1984), (2) they are generated in the ISM several kpc from a 
'naked' Seyfert nucleus (Korista & Ferland 1989). The second of these possibilities makes 
predictions regarding relative line luminosities (assuming solar metallicity). Specifically, this 
model predicts that [NcV] A3426 will have a comparable line width to [FcX], and be ~12 
times stronger, and makes further predictions for the ratio of the luminosities of the iron lines 
to the luminosity of the narrow component of H/3. While the widths of [FeX] and [NeV] A3426 
are comparable (2100±600km s"^ vs. 2500±800km s"^), [NeV] A3426 is approximately 40 
times stronger than [FeX] (we corrected the [FeX] 6375 flux from contamination from [01] 
6364 by assuming [01] A6364/[OI] A6300 = |). Furthermore, the L([FeX])/H/3, L([FeXI])/H/3 
and L( [FeXIV] )/H/9 ratios predicted by Korista & Ferland (1989) and those in our spectra 
all differ by factors of between 2 and 30, depending on the assumed Balmer decrement. The 
first of the models predicts that these iron lines may show a shift (blue or red) compared 
to the systematic redshift of the system, and that the FWHMs of the iron lines will be 
intermediate between the forbidden lines and the broad components of the permitted lines. 
Both these criteria are met, hence we conclude that the high excitation iron lines in F05189- 
2524 originate in a CLR between the narrow and broad line regions. 

Under the assumption that these high excitation iron lines are produced in photoionized 
gas (rather than shocks in radio jets), then the lines can be used to constrain the temperature 
and density of this gas within a few pc of the central ionizing source. The detection of [FeXIV] 
imphes (Ferguson, Korista & Ferland 1997) a range in Hydrogen column density of 5.75 <log 
n(H)(cm~^) < 8.0, and a peak electron temperature of ~ 160, OOOK. The detection of this 
line also implies that we have a line of sight to regions very close to the central ionizing 
source; if the hydrogen column density range and electron temperature are correct then 
[FeXIV] is expected to be produced only within ~ 0.2pc of the broad line region. 
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Analysis 

The complex lines of sight to the AGN intimated by the presence of the high ionization iron 
lines makes deriving reliable optical line diagnostics complex. The Hf3 line cannot be reliably 
fitted with more than a single Gaussian and single Lorentzian component. Attempts to fit the 
Ha line with more than a single component were similarly unsuccessful; when a broad and 
narrow component were simultaneously fitted then the flux of the broad component varied 
by ~ 30%, whereas the flux of the narrow component varied by up to an order of magnitude, 
depending on the values of the input parameters to specfit. We therefore assumed, solely for 
the purposes of deriving line diagnostics, [Nil] A6583/ifQ; (narrow) = 1.0. For the Hj3 line we 
used only the Gaussian component when measuring the [OIII] \5007 /Hf3. The Hf3 flux we 
used is thus likely too low as it does not include any contribution from a narrow Lorentzian 
component. The value of log([OIII] A5007/if/3) plotted in Figure 7 is thus probably too high, 
and indeed a value of 1.56 is higher than any other published AGN. This overestimate however 
should not affect the classification. The three primary line ratio diagnostics all classify this 
object as an AGN, as do the [OIII] X5007/HP vs. [Oil] A3727/[OIII] A5007 and [OIII] 
X5007/HP vs. [Oil] AA7320, 7330/ifQ; diagnostics. Furthermore, as the I band spectrum 
shows the [SIII] AA9069, 9531 hues we were able to employ the [SIII] AA9069, 9531/i/a vs. 
[Oil] AA7320, 7330/i/a diagnostic described in Osterbrock, Tran & Veilleux (1992), which 
also classifies this spectrum as an AGN. 

With an assumed intrinsic Balmer decrement of 3.10, the extinction corrected star for- 
mation rate derived using Equation 5 is 30 ± IOMq yr~^; this should of course be regarded 
as an upper limit due to the presence of an AGN. The [SII] A6716 & [SII] A6731 lines could 
not be reliably deblended in this object, hence we only present an estimate of the electron 
temperature in the narrow line region. For an assumed intrinsic Balmer decrement of 3.1 the 
[Oil] A3727 / [Oil] A7325 ratio is 16.98, which predicts an electron temperature of 13, 900K 
for an assumed electron density oi — 1 x lO^cm"^, and 6, 200K for an assumed electron 
density of ng = 1 x lO^cm"^. 

The strong, broad emission lines in the UV spectrum evidently arise from an AGN, 
however the analysis we could perform was limited, for two reasons. Firstly, most emission 
line diagnostics for the ionizing source in the narrow line region of AGN use one or both 
of the OVI A1034 and GUI] A1909 lines (Hamann & Ferland 1999), and both these lines he 
outside the G140L bandpass at 2; = 0.043. Secondly, due to the hmited spectral resolution of 
our data it proved impossible to separate the broad and narrow components of most of the 
emission lines in the UV; the only line this was achieved for was CIV A 1549 and even in this 
case the la errors on the resulting widths and fluxes were ~ 30%. We therefore first used 
the [OIII]A5007/H/? vs. [NeV]A3426/[NeIII]A3869 diagnostic (Allen, Dopita & Tsvetanov 
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1998) to distinguish between shocks and photoionization, which imphed a photoionizing 
continuum with a ~ — 1. Wc then used the NIV]/CIV diagnostic presented by Hamann et 
al (2002) and the total CIVA1549 flux, to derive a lower hmit to the metaUicity, which gave 
log(Z/ZQ) > 0.4. 

Overall, F05189-2524 knot 1 clearly contains an AGN, based on both the hne shapes in 
the optical and UV and the optical line diagnostics. Furthermore, we have a line of sight 
very close to the central regions of this AGN, based on the broad component on H/? and 
the detection of the high excitation 'coronal' iron lines. Although there are no unambiguous 
signs of a starburst in the optical spectra, the presence of NIV] A 1488 in the UV suggests a 
large population of W-R stars, implying that a young (3-4Myr) starburst is also present. 

Finally, we note that an interesting experiment would be to measure the polarization of 
the high-ionization lines found in this spectrum. The broad components found on the Ha 
and H/5 lines in our spectra, and results from previous authors (Young et al 1996; Veilleux, 
Sanders & Kim 1999) imply that the BLR in this knot is moderately, though not completely, 
obscured. As the distance of the source of the high-excitation iron hues from the BLR is 
predicted to decrease with increasing ionization stage (Ferguson, Korista & Ferland 1997), 
the relative polarizations of these lines would be an interesting diagnostic of the change in 
obscuration as a function of distance from the BLR. 



3.2.3. Knot two 

This knot is fainter than knot one, and is not detected in the UV. Due to the spatially 
extended emission lines in knot 1 , which contaminate most of the key diagnostic lines in knot 
2, little quantitative information can be derived. Nevertheless, a simple comparison of the 
optical spectra of knots 1 and 2 shows two clear differences. Knot 2 shows a distinct Balmer 
break and Na ID A5896 absorption, neither of which is seen in knot 1. Furthermore, knot 2 
does not show the Hell A4686 line; based on the H/3 and [Nelll] A3869 fluxes of the two knots 
we would have expected to detect Hell A4686 in knot 2, if it were present. We speculate that 
knot 2 has a higher obscuration level than knot 1, and somewhat lower ionization, though 
our data do not allow us to quantify this. Overall, the properties of this knot combined with 
the properties of knot 1 are not consistent with the central region of F05189-2524 being a 
single nucleus bisected by a dust lane, contrary to the results of Surace & Sanders (1999) 
and Scoville et al (2000). Instead, the most hkely interpretation is that F05189-2524 is a 
single nucleus system, and that knot 2 is a second line emitting region ionized by the nucleus 
in knot 1, but with a different level of excitation and obscuration. 
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One puzzling point is the clear detection of a line we identify as [FeX] 6375 in knot 2. 
Careful inspection of the 2-d image shows that this line does not arise via contamination 
from knot 1, but if it does arise from knot 2 then the apparent lack of other high excitation 
lines is strange. The simplest scenario which fits with a single nucleus interpretation is that 
the [FcX] emission from knot 2 is a result of scattering and reflection from knot 1, in which 
case we would expect this line to be polarized, but we cannot quantify this or explain the 
absence of the other high excitation iron lines without higher quality data. 



3.3. IRAS F12071-0444 

HST imaging shows a central, reddened, compact nucleus surrounded by a 'chaotic' 
high surface brightness region and a number of small blue knots, and two tidal tails (Suracc 
et al 1998). Ground based optical spectroscopy shows a Sy2 spectrum (Veilleux et al 1995; 
Kewley et al 2001b), and Veilleux, Sanders & Kim (1997) use near- infrared spectroscopy to 
infer the presence of a hidden BLR. 

A single spectrum is evident in the STIS B and 1 band frames, with nothing detected in 
the UV frame. The optical spectra contain several high excitation Neon lines and the Mgll 
AA2796,2803 doublet. Both H7 and RS are seen in emission. The optical lines are symmetric 
about the continuum, with spatial extents of 800 - lOOOpc (for F 12071-0444 the STIS pixel 
scale is ~ 120 pc pixel"^ in the optical). Two high excitation iron hues are also present 
in the I band spectrum, [Fe VII] A5721 and A6087 (Table 7), detected at the ~ 4(7 level. 
All the line ratio diagnostics classify this spectrum as an AGN, and as a Sy2 rather than 
a LINER. The [SII] AA6716, 6731 doublet was not resolved for this object, hence we only 
estimated electron temperatures in the narrow line region. The [Oil] A3727 / [Oil] A7325 
ratio, at 55.58, predicts a temperature in the NLR of 9700K for He = lOOcm"^, and 7200K 
for He — lOOOcm"^. Overall we classify this object as an AGN, with no clear signs of an 
accompanying star burst, though we cannot rule out the presence of a star burst in the other 
knots in this object. 



3.4. IRAS F15206+3342 

Sanders et al (1988b) classify this galaxy as a Sy2 on the basis of ground based optical 
spectroscopy. Veilleux, Sanders Sz Kim (1997) do not however find any evidence for a hidden 
BLR. HST and ground-based imaging shows a complex morphology, with a large number 
of knots and a tidal tail (Surace et al 1998; Veilleux, Kim, & Sanders 2002). Most of these 
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knots are probably starbursts, but two may be active nuclei. Arribas & Colina (2002) find 
evidence for gas inflow along the tidal tail, and postulate that this object is a late stage 
merger between two disk galaxies which suffered a retrograde encounter. 

For this object the STIS slit was oriented so that it covered the centers of two knots, 
one of which was identified as a possible AGN by Surace et al (1998), and partially covered 
a third (Figure 1). Three spectra are thus visible in the STIS frames, and are discussed in 
the following subsections. 

3.4-1- Knot one 

Emission & absorption lines 

Knot 1 is the brightest of the three in the UV, and the UV spectrum contains a number of 
absorption systems, including Lja. Only two emission lines are present; NV A1240 and CIII 
A1175. The (unsaturated) absorption systems have FWHMs of 800 - 1500 km s~^ and show 
a blueshift of ~ 300km s~^. In addition to the commonly seen UV absorption hues there 
are two weak features at observed-frame 1598A and 1609A that proved difficult to identify; 
features at these wavelengths are (to our knowledge) not seen in spectra of AGN or starburst 
galaxies, or composite QSO spectra, and do not correspond to absorption features commonly 
seen in Damped Lyman Alpha systems (Prochaska et al 2001). These two features arc most 
likely spurious, and so are not marked in Figure 6, even though they are detected at ~ 3.5(T. 
If they are real, then the most likely candidates, based on currently available atomic line list 
data, arc multiplets of allowed and semi-forbidden transitions of Fe II and Fe III (another 
possible ID motivated by absorption features seen in DLAs for the longer wavelength feature 
is Pbll A1433), but it is beyond the scope of this paper to provide a secure identification. 

The optical spectra show a steep continuum and a Balmer break. The emission lines, 
which are offset from the continuum by ~ 300pc, have a spatial extent of 800 - 900pc. We 
detect both the [NeV] A3426 and [Nelll] A3869 lines, but no high excitation iron hues. A 
number of absorption features are present blueward of 4000A; we identify Mgll AA2796,2803 
and (cautiously) Mgl A2852, Hel A2945, Fell A3180 and Fel+Fell A3261. These absorption 
features show blueshifts of between 300km s"^ (for Mgll AA2796,2803) and 2600km s"^ (for 
Fel+Fell A3261) relative to the systemic redshift as measured from [OII]A3727. The optical 
emission hue widths are narrow, with a corrected [OIII]A5007 FWHM of 620 ± 50 km s~^. 
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Analysis 

It proved impossible to disentangle the Ha and [Nil] emission; attempts to do so produced 
errors of at least a factor of two on the derived Ha flux depending on the input paramaters 
to specfit, hence the same procedure as for IRAS F01003-2238 was used to derive line ratio 
diagnostics. There is some ambiguity in the classification of this spectrum. The [OIII] 
A5007////3 vs. [Nil] A6583/i/Q; and [OIII] A5007/i//3 vs. [SII] AA6716, 6731///q; diagnostics 
classify this knot as either an AGN or a star burst, depending on the assumed intrinsic 
Balmer decrement. The weak upper limit on [01] A6300 renders diagnostics based on this 
line inconclusive. The [OIII] X5007/H(3 vs. [OH] AA7320, 7330/i/a and [OIII] \5007/H(3 
vs. [OH] A3727/[OIII] A5007 diagnostics classify this object as an AGN and HII region 
respectively, irrespective of any assumptions. We do not present estimates for electron 
densities or temperatures in the NLR as the [OH] A7325 line is not detected, and the [SII] 
doublet could not be deblended. With an assumed intrinsic Balmer decrement of 2.85, the 
extinction corrected star formation rate (if all the [Oil] A3727 flux arises from star formation) 
derived using Equation 5 is 60 ± ISM© yr"-*^. 

To quantify the presence of a starburst and/or AGN within this knot, we have modeled 
the UV spectrum using the Starhurst99 v4.0 code (Leitherer et al (1999) and the same 
parameter ranges as for F01003-2238 (§3.1.4), under the ah initio assumption that the UV 
spectrum is starburst dominated. The observed UV spectrum is not consistent with a single 
episode or any combination of two episodes of star formation (Figure 12), irrespective of 
the adopted age, IMF slope and lower stellar mass limit, or whether the star formation is 
continuous or a 'burst'. The emission line we identify as CHI 1175A is not present (this 
line is in fact not predicted to exist in any of the starburst models we investigated), and 
the multiple strong absorption features in the observed spectrum are poorly matched by the 
model. Significant absorption at 1465A due to 01 A1304 and Sill A1307 could be generated 
by a starburst with a very steep IMF slope (> 7) but then the other absorption lines, 
particularly Lya, would not be consistent with the data. A close match around the Lya and 
NV fines can however be obtained with a single starburst model. The remainder of the UV 
spectrum requires a different explanation. 

There are two possible origins for the absorption features in the UV spectrum of 
F15206+3342. The first is that they arise due to interstellar absorption. Interstellar ab- 
sorption features in the rest-frame UV that are qualitatively similar to those seen in our 
spectrum are seen in, for example, UV spectra of Lyman Break Galaxies (Pettini et al 2000; 
Shapley et al 2003). It is possible therefore that the features in our spectrum are due to 
interstellar absorption. An argument against this however is that most of the absorption 
features in our spectrum are between 2 and 5 times broader than those seen in Lyman Break 
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Galaxy spectra. For example, for the SilV AA1394, 1403 doublet; Pettini et al find rest-frame 
widths of I.94A and 1.42A, Shapley et al find widths of l.SSA and O.8IA, whereas we find 
widths of 4.65A and 3. 71 A. The exception is the CII A1335 line, which has a comparable 
width in our spectrum and that of Shapley et al. 

The second possible origin for these features is that they arise all or in part from a Broad 
Absorption Line (BAL) region, itself an outflow driven by a hidden AGN. Approximately 10% 
of QSOs show BALs with widths > 2000km (Weymann et al 1991; Arav ct al 2001; Hall et 
al 2002; Lacy et al 2002; Reichard et al 2003). The existence of an AGN driven BAL outflow 
in F15206+3342 is supported by the presence of Mgll AA2796,2803 and other absorption lines 
seen in BALs (Hall et al 2002) in the B band spectrum, which would classify this object as a 
low ionization BAL, or 'LoBAL', which comprise ~ 10% of the BAL population. The UV and 
B band spectra contain absorption features that are thought to arise from metastable excited 
levels of iron; such features are very rare amongst LoBALs, and would place F15206-I-3342 
knot one in the so-called 'FeLoBAL' class of AGN (Hazard et al 1987; Becker et al 1997; 
Lacy et al 2002). This picture however is also problematic, for two reasons. Firstly, the 
absorption lines are narrow compared to the classical deflnition of absorption lines in a BAL 
QSO (Weymann et al 1991); the corrected FWHMs of Mgll AA2796,2803 and SilV A1394 
in our spectra are 1500 ± 150km s~^ and 1100 ± 200km s~^, though Hall et al (2002) do 
find absorption troughs of comparable widths to ours in some SDSS BAL QSOs. Our object 
would thus be formally termed a 'mini- FeLoBAL' AGN (Hamann 2001). Secondly, a BAL 
nucleus cannot be responsible for all the observed features in the UV spectrum. BAL QSOs 
show absorption in GIV, NV, SilV and Lya in order of decreasing strength (Hall et al (2002) 
& rcfs therein) whereas our spectrum shows strong Lya absorption and weak NV emission. 
Overall therefore, we conclude that the most likely scenario is that the Lya absorption and 
NV emission arise from a starburst, and the UV and optical absorption features arise from 
an outfiow driven by a mini-FcloBAL AGN, but we cannot discount the possibility that some 
or all of the absorption features (particularly Gil A1335) arise due to interstellar absorption 
of a starburst wind. 

If the BAL AGN interpretation is correct then it would have interesting implications 
for the origin of the BAL QSO phenomenon. There are currently two models for the origin 
of the absorption lines in BAL QSOs; the 'disk wind' model and the 'youth' model. In the 
disk wind model (Murray & Ghiang 1998; Proga, Stone & Kallman 2000) a wind from an 
accretion disk encounters a high column density, highly ionized gas. The high ionization 
state of this gas means it can transmit UV photons from (and thus be driven outwards 
by) resonance lines such as GIV A 1549, but the high column density of the gas shields it 
from higher energy photons that would otherwise completely ionize the gas. In this model 
BAL QSOs are those QSOs viewed along a line of sight that coincides with this outfiowing 
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gas, marking the detection of BALs in QSOs as largely orientation dependent (Elvis 2000). 
The 'youth' model on the other hand contends that BAL QSOs arc young objects, still 
surrounded by large quantities of gas and dust; in this case the BALs do not arise due to a 
particular line of sight (Voit, Weymann & Korista 1993; Becker et al 1997; Williams, Baker 
& Perry 1999). The 'youth' model is often applied to LoBAL QSOs; Sprayberry & Foltz 
(1992) propose that the differences between the line properties of LoBAL QSOs and those 
of non-BAL QSOs cannot be explained solely in terms of different relative orientations. 
Voit, Weymann & Korista (1993) propose an alternative scenario in which LoBALs form 
via ablation of dust by UV photons in outflows arranged either as a thick disk or as an 
isotropic distribution of clouds; in this case LoBALs do fit within AGN orientation schemes, 
and this scenario is consistent with polarimetric observations of large samples of BAL QSOs 
(Schmidt & Hines 1999; Ogle et al 1999; Hutsemekers & Lamy 2000). Even for this model, 
Voit, Weymann & Korista (1993) note that BAL QSOs would likely evolve into less obscured 
QSOs as the outflowing material either settles into a thinner conflguration or entirely ablates 
away, and propose that this scenario is consistent with LoBAL QSOs being youthful objects. 
An evolutionary hnk between LoBAL QSOs and ULIRGs has also been previously suggested 
from other observations (Lipari, Colina & Macchetto 1994; Canahzo & Stockton 2001). The 
presence of a BAL nucleus within F15206+3342 would add further credence to the 'youth' 
model for LoBAL QSOs. 



34.2. Knot two 

Knot 2 is the faintest of the three in the optical, and the UV emission is also very weak. 
This is due at least in part to the positioning of the sht; as can be seen from Figure 1 the 
slit passes over only one edge of this knot. The discussion of this spectrum should therefore 
be regarded with caution. In the optical, similar emission lines are present as in knot 1, 
including the [Nelll] A3869 line, but also the H7 and Hel A5876 lines. The [NeV] A3426 
line may be double peaked; adopting the systemic redshift then the rest-frame wavelengths 
of both features, at 3396A and 3422A lie just within the range seen in (for example) SDSS 
QSOs (Vanden Berk et al 2001), but we conservatively only identify the longer wavelength 
feature as [NeV] A3426. No stellar absorption features were reliably identified. 

It proved impossible to disentangle the Ha and [Nil] emission, hence the same procedure 
was used as for knot 1. The results from our three primary diagnostics arc ambiguous between 
an HII and AGN, as is the [OIII] A5007/i7/3 vs. [Oil] A3727/[OIII] A5007 diagnostic, though 
the [OIII] \hm/Hl3 vs. [Oil] A3727/[OIII] A5007 diagnostic robustly classifies this object 
as an HII region. Overall therefore, we classify this knot as a starburst. 
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The [Oil] A7325 line is not detected in this object, however the [SII]AA6716, 6725 doublet 

was successfully deblended, hence we present estimates for the electron density in the NLR. 
The [SII] A6716 / [SII] A6731 ratio, at 0.90, infers electron densities of 1200cm-3, QlOcm-^, 
and 690cm"^ for temperatures of 5,000K, 10,000K and 20,000K respectively, which implies 
a moderately dense NLR. With an assumed intrinsic Balmer decrement of 2.85, the optical 
star formation rate derived using Equation 5 is 75 ± 2OM0 yr~^ 



Knot 3 is optically the brightest of the three, but is not detected at all in the UV. The 
lines arc symmetric, with a spatial extent of ~ 850pc. The line widths are narrow compared 
to the other spectra in the sample, with a corrected [0111] A5007 FWHM of 450 ±50 km s^^. 
Though the continuum detection was strong in both the B and I band spectra, no stellar 
absorption features were identified. 

For this object the lines are sufficiently narrow that the Ha and [Nil] lines can be reliably 
deblended. All five of our adopted line diagnostics classify this spectrum as a starburst rather 
than an AGN. Both the [Oil] AA3727, 7325 doublets arc detected, and the [SII]AA6716, 6725 
doublet was successfully deblended. Prom the [Oil] A3727 / [Oil] A7325 ratio, 31.53, we 
derive electron temperatures of 13900K, 9300K and 5100K for densities of lO^cm"^, lO^cm"^ 
and lO^'cm-^^ respectively. Prom the [SII] A6716 / [SII] A6731 ratio of 0.77 we derive electron 
densities of 1200cm~^, 1600cm~^ and 2100cm~^ for assumed temperatures of 5000K, lOOOOK 
and 20000K respectively. 

As this knot is universally classified as an HII region, and because the [Nil] and Ha 
fines could be deblended, we have estimated the oxygen abundance using the [Nil] A6584 
/[Oil] AA3726,3729 diagnostic presented by Kewley & Dopita (2002), though we note that 
this diagnostic is only reliable for metallicities above one half solar. We derive: 



34.3. Knot three 




(6) 



which corresponds to an oxygen abundance of: 



log{0/H) + 12 = 8.91 



(7) 



This oxygen abundance is higher than solar (8.69, AUende Prieto, Barklem & Asplund 
(2001)) and is comparable to the most metal rich local UV selected galaxies (e.g. Contini et 
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al (2002)). 

With an assumed intrinsic Balmer decrement of 2.85, the extinction corrected star for- 
mation rate derived using Equation 5 is 1000 ± SOOMq yr~^. Even by the standards of local 
ULIRGs, this star formation rate is extraordinarily high (Farrah et al 2003), and is compa- 
rable to that seen in distant HLIRGs (Farrah et al 2002b). If we instead use the extinction 
corrected Ha luminosity, and the calibration given by Kennicutt (1998): 

SFR = 7.9 X lO-^^L(Ha) (8) 

then we obtain a comparable value; 800 ± 25OM0 yr^^. Even if all the star formation in 
Fl 5206+3342 was in this knot and visible in the optical, this star formation rate is likely to 
be too high by a few hundred Mq yr^^. Establishing whether or not this is the case, and the 
reasons behind any overestimate, is beyond the scope of this paper, but we briefly mention 
two possibilities. As noted by Kennicutt (1998), star formation rate calibrations for Ha and 
[OII]A3727 depend strongly on both the optical extinction, and the assumed stellar IMF. It 
is quite possible that the measured extinction in this knot is incorrect. It is also possible 
that the IMF slope is steeper than that assumed in deriving Equations 5 and 8, this would 
mean there are a larger number of high-mass stars present, producing a greater ionising flux 
per Mq. 

The unambiguous classification of this knot as an HII region means that an AGN is 
unlikely to contribute significantly to the continuum emission. We can therefore estimate 
the stellar mass of this knot by combining the extinction measure from our observations with 
diagnostics from HST WFPC2 imaging (Suracc et al 1998); we note this is the only knot 
in our sample for which this analysis is possible, as the others arc cither weakly detected or 
contain an AGN. This gives a mass of ~ 2 x lO^M© (note this differs from the mass presented 
in Table 2 as we have applied a reddening correction based on our spectra). This translates, 
assuming the knot is approximately spherical, to a stellar mass density of ~ 2OOM0pc~^. 
Both the mass and stellar density of this knot are comparable to, but at the lower ends of, the 
ranges seen in the cores of elliptical galaxies, which typically have masses up to ~ lO^'^M© 
and stellar densities in the range 10^ — IO^Mq (Lauer 1989). Coupled with the very high star 
formation rate (at least a few hundred Mq yr~^) this makes this knot an excellent candidate 
for the forming core of an elliptical galaxy. 
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4. Sumrriciry Conclusions 

We have presented HST STIS longslit optical and UV spectra of the knots in the nuclear 
regions of four 'warm' ULIRGs, and used a combination of line ratio diagnostics and spectral 
synthesis modelling to constrain the active power source in each knot. A summary of the 
properties of the knots in each object is as follows: 

1) IRAS F01003-2238 - The STIS data show a single bright narrow line spectrum. The UV 
spectrum contains signatures of O supergiants, both WC and WN type Wolf-Rayet stars, 
and (probably) an AGN. Results from spectral synthesis modeUing are consistent with a 
'burst' of star formation with age 3-4Myr and an IMF slope < 3.3. The optical spectra 
show what may be high excitation iron lines, indicating a line of sight close to the broad line 
region, and evidence for a wind, probably starburst driven, that originated in the nucleus 
around 0.5Myr ago. 

2) IRAS F05189-2524 - Two spectra can be seen, corresponding to the two sections of the 
central bisected 'nucleus'. The UV spectrum of knot 1 shows several strong permitted hnes 
from an AGN, together with weak NIV] A1488 emission that probably arises from a young 
(~ 4Myr) 'burst' of star formation. The B and I band spectra of knot 1 show a large 
number of lines, including the OIII A3 133 Bowen resonance- fluorescence line. Also present 
are several high excitation iron lines, including [FeXIV] A5303, which imply lines of sight 
from the narrow line region through to ~ 0.2pc from the accretion disk. The optical line 
diagnostics classify this knot as an AGN. Further diagnostics show no evidence for shock 
heating, and imply a metallicity of 2.5Z0 or higher. Knot 2 is fainter than knot 1, and is 
not detected in the UV. Its properties are consistent with a higher obscuration and lower 
ionization level than knot 1. Considered together, the properties of the two knots are not 
consistent with a single nucleus bisected by a dust lane, but rather with a single nucleus 
(knot 1) together with a second knot of line emitting gas ionized by the nucleus in knot 1 
but with different ionisation conditions and obscuration. 

3) IRAS F12071-0444 - A single spectrum of the nucleus can be seen. This object is unde- 
tected in the UV, and the continuum is only weakly detected in the B and I bands. Only 
a few optical emission lines are present, limiting the scope of the analysis. The optical line 
ratio diagnostics universally classify this object as a Sy2. There is no evidence for an ac- 
companying starburst, though we cannot rule out the possiblility of starbursts in the other 
knots in this object. 

4) IRAS F15206+3342- Three spectra can be seen. The UV spectrum of knot 1 is a chimera; 
containing Lja and NV lines that can be modeled as arising from a starburst, together with 
a CIII A1175 line that is not present in any spectral synthesis model we investigated, and 
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thus likely arises from an AGN. The UV and B band spectra contain a number of absorption 
systems that probably arise from a hidden FeLoBAL AGN, the alternative being that they 
arise from interstellar absorption. Knot 2 is only weakly detected; this spectrum is classified 
as an HII region with a moderately dense narrow line region. Knot 3 is optically bright 
but not detected in the UV. Line ratio diagnostics classify this knot as an Hll region with 
a supersolar (Z^l.SZ©) metallicity and a very high star formation rate. Combined with 
previous results this implies that this knot is plausibly the forming core of an elliptical 
galaxy. This is the only object in the sample for which the knot classifications given by 
Surace et al (1998) were found to be incorrect, for knots 1 and 3. 

The knots in our sample of 4 'warm' ULlRGs are remarkable both for their similarities 
and their differences. Three of the four objects contain a knot with clear AGN signatures, 
and one object (F01003-2238) also hkely contains an AGN, demonstrating that selection 
based on 'warm' mid-IR colors does preferentially find ULIRGs with an AGN. In three out 
of four cases this AGN was found in the knot identified as a putative AGN from multiband 
HST imaging. Three of the four ULIRGs in our sample also contain at least one starburst 
knot, confirming that most ULIRGs are 'composite' objects, containing both a starburst 
and an AGN. The detection of luminous AGN and/or starbursts in all of the knots in our 
sample constitutes plausible, though not compelling evidence that these knots are the sites 
for the dust-shrouded starburst and AGN activity that power the IR emission in ULIRGs. 
The AGN display a wide range of properties; two are classical narrow line AGN, one shows 
both broad and narrow lines and has multiple lines of sight to the center, and one is a narrow 
line AGN with (plausibly) an FeLoBAL outflow. Of the starbursts, only 3 have constraints 
on their ages, with two having ages 3-4Myr and the third with an age of < 20Myr. Though 
this could be a selection effect, it seems that young bursts of star formation are common in 
ULIRGs, indicating that most ULIRGs probably undergo several such bursts during their 
lifetime. Finally, the properties of one starburst knot arc consistent with it being the forming 
core of an elliptical galaxy, the first time to our knowledge that such a structure has been 
identified in a ULIRG. The identification of such a structure is supportive of the idea that 
ULIRGs can evolve into systems similar to local elliptical galaxies. 

Our study has demonstrated the power of combining high spatial resolution optical and 
UV spectroscopy to study local ULIRGs; the optical data, in addition to providing fines 
for standard line ratio diagnostics, also isolated several rare, high excitation lines in most 
systems. In several cases the AGN was only unambiguously seen in the UV, and quantitative 
constraints on the properties of the starbursts only came from modelling the UV spectra. 
In one case (F01003-2238) this approach showed (probable) AGN signatures, when such 
signatures had previously only been seen in mid-IR spectroscopy. While this method does 
not directly measure what contributes to the IR emission it does give a sensitive, largely 
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model free measure of what co?//c? contribute to the IR emission. Indeed, this study would not 
have been possible without a space-based optical/UV spectrograph with spatial resolution 
of < 0.1". Future multiwavelength studies focused on these knots are likely to give powerful 
insights into the nature of starburst and AGN activity in ULIRGs, and into ULIRG evolution. 
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Table 1. 'Warm' ultraluminous infrared galaxy sample 



Galaxy 


RA (J2000) 




Dec 


Redshift 


r . a 


ruB 


IRAS F01003-2238 


01 


02 


49.99 


-22 


21 57.5 


0.1177 


12.39 


18.5 


IRAS F05189-2524 


05 


21 


01.47 


-25 


21 45.4 


0.0426 


12.20 


15.6 


IRAS F12071-0444 


12 


09 


45.12 


-05 


01 13.9 


0.1284 


12.51 


18.2 


IRAS F15206+3342 


15 


22 


38.05 


+33 


31 35.9 


0.1244 


12.33 


16.8 



"Logarithm of the 8 — 1000//m luminosity, taken from Kim & Sanders (1998) and rescaled 

to our cosmology. 



Table 2. Knot magnitudes, radii and masses 



Galaxy 


Knot 


rriB 


mi 


Radius 


Mass 










pc 




F01003-2238 




19.27 


18.32 


28 


3 X 10* 


F05189-2524 




18.78 


16.55 


25 


2 X 10^° 




2 


19.01 


17.37 


27 


1 X 10^ 


F12071-0444 




21.59 


19.00 


244 


4 X 10^0 


F15206-I-3342 


1 


19.45 


18.54 


88 


3 X 10^ 




2 


19.50 


18.68 


104 


2 X lO'^ 




3^ 


19.29 


17.37 


57 


3 X 10^°^ 



Note. — All data are taken from Surace et al (1998). 
Knots are numbered as in Figure 1. 

^knots identified by Surace et al (1998) as putative AGN 
based on their broad-band colors. 

'^stellar mass inferred by Surace et al (1998) assuming 
that all of the B and I band flux is due to stellar emission 



'^but see §3.4.3 
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Table 3. Ultraviolet emission line fluxes 



T * 

Line 


Wavelength 

A 


F01003-2238 
1 


F05 189-2524 
1 


F12071-0444 
1 


F15206-I-3342 
1 


rieii 


iUoO 


1 1 Q 

i.iy 








cm 


1175 


abs 


0.13: 




1.23: 


Lya 


1216 


10.83 


52.74 




abs 


NV 


1240 


1.58: 


20.51 




0.22: 


Sill 


1260 


abs 






abs 


CII 


1335 




0.90 




abs 


SilV 


1394,1403 


1.08: 


2.74 




abs 


NIV] 


1488 


0.53: 


1.12 






CIV 


1548,1551 




28.54 






[NeV] ^ 


1575 




0.30:: 






[NelV]^ 


1601 




0.35:: 






Hell 


1640 




13.22 






Note. 


— Fluxes are 


in units of 10" 


-15gj,grs Cm~^S" 


and have not been corrected 



for extinction. Flux errors are of the order 5% — 10% for most sources; those fluxes 
marked with a colon have flux errors of ~ 20% and those fluxes with a double colon 
have errors of ~ 30%. 

'^Uncertain ID 
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Table 4. Optical emission line fluxes 



Line 


Wavelength 
A 


F01003-2238 
1 


F05189-2524 
1 2^ 


F12071-0444 
1 


F15206+3342 
1 2 3 


iVigii 


07QR OSIT^ 


O.lO 


N/A 


N/A 


n 49- 

U.4Z. 


abs 






A/TrrT 

Mgi 


ZoOZ 




N/A 


N/A 




abs 






TTbTT 

reii 


zyoo 




N/A 


N/A 




abs 






nei 


OQ/I ^ 

zy^to 




N/A 


N/A 




abs 






DTTT 


100 




7.13 


N/A 










reii 


OloU 






N/A 




abs 






TToT 
nei 


oioy 




5.26: 


N/A 










T?r.Tj_T?QTT 
rCi-hreii 


OZOl 






N/A 




abs 






[i N e V J 


004D 




21.96 


N/A 










[In e V J 


o4ZD 


Z.Do. 


58.44 


0.72: 


u.oy. 


0.41: 


0.34: 


0.55 


[UllJ 


/ z / 


1 97- 
l.Z / . 


3.91: 




n 71 
u. / 1 


1.70 


1.55 


2.01 


\L\ eiiij 


oooy 


fi no 
o.uy 


17.11 


0.38: 


1.10 


0.52:: 


0.16:: 


1.14 


riMcTTTl 

[iNeiiij 


oyoo 


9 9R 
Z.ZD 


5.50 


abs 


n Q7- 








[bllj 


4U0o,4U ( D 




1.75: 












no 


41UZ 


0. io.. 


2.09 


abs 


n /1 1 •• 






0.63: 


117 






5.51 




n ztn- 




0.43 


1.84 


1\TTTT 
IN ill 


/1R/in AP.AO 
4D04,4U4U,4D4Z 


l.ZD. 














TTpTT 


4L)OL) 


x.zo.. 


1 Q9 
10. yz 




U.L)4 












4 Q6 


34.12 


0.23: 


1 49 


0.95 


1.25 


5.49 


[UiiiJ 


4yoy 


Q 07 

y.u / 


62.44 


N/A 


D.yi 


1.28 


1.32 


8.37 


[OIII] 


5007 


23.67 


177.98 


N/A 


20.82 


4.20 


4.31 


21.43 


Hell 


5412 




5.72 












He! 


5876 


_b 


4.28 








0.35: 


2.66: 


[01] 


6300 


3.23 


4.37 




1.97 








[Nil] 


6548 




13.30 


N/A 


6.17 






6.54 


Ha 


6563 




202.20 


N/A 


13.99 






68.26 


[Nil] 


6583 




39.94 


N/A 


18.54 






19.66 


Hq;+[NII] 




56.94 




N/A 




8.76 


12.49 




[SII] 


6716,6731 


3.29 


12.95 


0.51: 


5.94 


1.71: 


1.15 


3.76 


Hel 


7065 




1.47: 
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Fig. 1. — HST WFPC2 images of our sample, taken from Surace et al (1998), with the sht 
positions superimposed and the nuclear regions marked. North is up and east is to the left. 
The scale bar denotes IKpc. 



Table 4 — Continued 



Line 


Wavelength 


F01003-2238 


F05189-2524 


F12071-0444 


F15206+3342 




A 


1 


1 


2* 


1 


12 3 


[Arlll] 


7138 




5.14 




1.29 


- - 1.96:: 


[Oil] 


7320,7330 


2.28: 


2.86: 




0.58: 


- - 2.27: 


[SIII] 


9069 




17.22 


N/A 




- - 8.35:: 


[SIII] 


9531 




25.08: 


N/A 







Note. — Fluxes arc in units of 10~^^ergs cm~^s~^ and have not been corrected for 
extinction. '-': Undetected. 'N/A': line is either outside the bandpass, contaminated 
by companion spectra, or in a noisy part of the spectrum. Flux errors are of the 
order 5% — 10% for most sources; those fluxes marked with a colon have flux errors 
of ~ 20% and those fluxes with a double colon have errors of ~ 30%. 

^extracted using a 1 pixel (42pc) aperture, fluxes are uncertain due to contamina- 
tion from knot 1. 

•^but see notes for Table 7 
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Table 5. Corrected Full Widths at Half Maximum of selected emission lines 



Line 


F01003-2238 


F05189-2524 


F12071-0444 


F15206+3342 




1 


1 


1 


1 


2 3 


Lya 


2200 ± 130 


N/A 




abs 




CIV 1549 




4600 ± 1000(B) 
1000 ± 500 (N) 








[Oil] 3727 


860 ± 500 


900 ± 275 


500 ± 100 


600 ± 300 


600 ± 50 550 ± 100 


H/3 


2100 ± 170 


3700 ± 900 (B) 
800 ± 100 (N) 


525 ± 100 


600 ± 100 


500 ± 50 450 ± 50 


[OIII] 5007 


2400 ± 70 


2000 ± 300 (B) 
900 ± 50 (N) 


950 ± 50 


650 ± 50 


550 ± 50 450 ± 50 


[FeXIV] 5303 




4800 ± 500 








[Fe VII] 5721 


900 ± 500*^ 


3000 ± 200 


900 ± 700 






[Fe VII] 6087 


900 ± 500^ 


2200 ± 100 


400 ± 700 






[01] 6300 


2200 ± 400 


950 ± 300 


1100 ± 150 




600 ± 300 


[Fe X] 6375 




2000 ± 600 


800 ± 350 







Note. — Widths are given in km s~^. FWHM of Lyct and NV cannot be measured in F05189- 
2524 as the broad and narrow components of these two lines cannot be rehably deblended, see text 
for details. 

*but see notes for Table 7 
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Table 6. Optical line ratios and spectral classifications 



Galaxy 




Ha 


Ay 


[OIII] 


[077]3726 [Nil] 
[OIII] Ha 


[SII] 

±1 Ot 


[OI] 

±1 OL 


[077]7325 
Ha 


(1) 


(2) 


(3) 


F01003-2238 


1 


8.61 


3.46 


4.77 


0.05 


0.25 


0.08 


0.05 


0.05 




- 


- 


- 






2.85 


0.0 


4.15 


0.20 


0.25 


0.07 


0.06 


0.04 




HII 


HII 


AGN 


F01003-2238 


1 


4.92 


1.45 


4.77 


0.05 


1.0 


0.13 


0.09 


0.09 




- 


- 


- 






3.10 


0.0 


4.50 


0.09 


1.0 


0.13 


0.10 


0.08 




AGN 


HII 


AGN 


F05189-2524 


1 


35.08 


7.87 


39.08 


0.02 


0.25 


0.08 


0.03 


0.02 




- 


- 


- 






2.85 


0.0 


28.42 


0.46 


0.25 


0.07 


0.04 


0.01 




AGN 


AGN 


AGN 


F05189-2524 


1 


8.77 


3.26 


39.08 


0.02 


1.0 


0.32 


0.11 


0.07 




- 


- 


- 






3.10 


0.0 


34.25 


0.08 


1.0 


0.30 


0.12 


0.05 




AGN 


AGN 


AGN 


F12071-0444 


1 


9.87 


3.63 


14.69 




1 .OO 


0.42 


0.14 


0.04 




- 


- 


- 






3.10 


0.0 


12.68 


0.38 


1.31 


0.39 


0.16 


0.03 




AGN 


AGN 


AGN 


r i5zUD+oo4z 


i 


0.88 


Z.76 


A /in 
4.40 


0.40 


0.25 


0.z6 


<0.18 


<0.11 












2.85 


0.0 


3.93 


1.17 


0.25 


0.25 


<0.09 


<0.0 


B 


HII 


A 


A 




1 


3.93 


0.75 


4.40 


0.40 


1.0 


0.46 


<0.20 


<0.20 












3.10 


0.0 


4.27 


0.54 


1.0 


0.45 


<0.20 


<0.18 


AGN 


AGN 


A 




2 


7.50 


3.03 


3.45 


0.36 


0.25 


0.12 


<0.03 


<0.03 












2.85 


0.0 


3.05 


1.16 


0.25 


0.11 


<0.04 


<0.02 


HII 


HII 


HII 




2 


4.29 


1.02 


3.45 


0.36 


1.0 


0.21 


<0.06 


<0.06 












3.10 


0.0 


3.31 


0.53 


1.0 


0.21 


<0.06 


<0.05 


AGN 


HII 


A 




3 


12.44 


4.62 


4.40 


0.08 


0.29 


0.06 


<0.05 


0.03 














2.85 


0.0 


3.65 


0.49 


0.28 


0.05 


<0.02 


0.02 




HII 


HII 


HII 



Note. — (1) [OIII] A5007/if/5 vs. [Nil] \6583/ Ha classification. (2) [OIII] \5007/H(3 vs. 
[SII] AA6716,6731/ifa classification. (3) [OIII] X5007/ H(3 vs. [01] A6300/iya classification. 'HII': 
Hll/Starburst spectrum. 'AGN': AGN spectrum. 'A': ambiguous classification between an HII and AGN 
spectrum. For F01003-2238 and F15206+3342 knots 1 & 2 the [Nil] and Ha lines could not be deblended, 
and for F05 189-2524 the narrow and broad components on Ha could not be deblended, hence two sets of 
line ratio diagnostics have been derived for these three objects assuming intrinsic Balmer decrements of 
2.85 and 3.10 respectively. In these cases the ratio is shown in bold. Uncertainties on the line ratios 
are ~ 10%, except for those ratios marked with a colon, where the uncertainty is ~ 20%. 
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Table 7. High excitation iron lines 



Line 


wavelength 
A 


F01003-2238 
1 


F05189-2524 
1 2^ 


F12071-0444 
1 


[Fe VI] 


5530 




2.07: 




[Fc VI] 


5632 


1.95^ 


3.66 




[Fe VII] 


3582 




3.32 




[Fe VII] 


3760 




6.39: 




[Fe VII] 


5159 




16.15 




[Fe VII] 


5721 


1.5:'^ 


19.96 


0.46: 


[Fe VII] 


6087 


1.00:*^ 


31.76 


0.58: 


[Fe X]^ 


6375 




9.44 0.53 




[Fe XI] 


7892 




7.85 




[Fe XIV] 


5303 




18.58 





Note. — ine fluxes are in units of 10~^^ergs cm~^s~^ and have not 
been corrected for extinction. Flux errors are of the order 5% — 10% 
for most sources; those fluxes marked with a colon have flux errors of 
~ 20%. 

^extracted using a 1 pixel (~ 42pc) aperture, IDs and fluxes are un- 
certain due to contamination from knot 1, see §3.2.3 for more details. 

^corrected for contamination from [01] A6364 assuming [01] 
A6364/[OI] A6300 = |. 

^Uncertain ID, may be blended with CIIIA5696, [NII]A5756 

^Uncertain ID, may be blended with CIVA5808, HeIA5876 

^Uncertain ID 
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Fig. 2.— G140L, G430L and G750L spectra of IRAS F01003-2238. The spectra have been 
smoothed with a 3 pixel boxcar and are plotted in the observed frame. Line IDs for the iron 
features at (observed-frame) 6416A, 6535A and 694lA are uncertain, see §3.1.2 for details. 
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Fig. 3.— G140L, G430L and G750L spectra of IRAS F05189-2524 knot 1. The spectra have 
been smoothed with a 3 pixel boxcar and are plotted in the observed frame. 
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Fig. 4.— G430L and G750L spectra of IRAS F05189-2524 knot 2. The spectra were extracted 
using a 1 pixel aperture to minimize contamination from the brighter knot, and only those 
lines with no contamination from knot 1 are marked. The spectra have been smoothed with 
a 3 pixel boxcar and are plotted in the observed frame. 
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[Olll] 



Mgll 



[Oil] 

[NeV] I [Nell 
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III] / 
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[Olll] 
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Wavelength (nm) 
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[Arlll] 



700 aoo 
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Fig. 5.— G430L and G750L spectra spectra of IRAS F12071-0444. The spectra have been 
smoothed with a 3 pixel boxcar and are plotted in the observed frame. 
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Fig. 6.— G140L, G430L and G750L spectra of IRAS F15206+3342 knot 1. The spectra have 
been smoothed with a 3 pixel boxcar and are plotted in the observed frame. The narrow 
'absorption' features in the I band spectrum are cosmic ray residuals. 
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Fig. 7.— G430L and G750L spectra of IRAS F15206+3342 knot 2. The spectra have been 
smoothed with a 3 pixel boxcar and are plotted in the observed frame. 
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Fig. 8.— G430L and G750L spectra of IRAS F15206+3342 knot 3. The spectra have been 
smoothed with a 3 pixel boxcar and are plotted in the observed frame. 
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log(NII/Ha) 




-1.2 -1.0 -0.8 -0.6 -0.4 

los(SII/Ha) 




Fig. 9. — Line ratio diagnostic diagrams. All fluxes used in these ratios have been corrected 
for the effects of extinction using the methods described in §2. Filled symbols show ratios 
where all lines have measured values, outline symbols connected by a line represent those 
objects where Ha and [Nil] could not be resolved and where two intrinsic values for the 
Balmer decrement (3.1 and 2.85) were therefore assumed. Arrows indicate upper limits. 
The curved solid lines are taken from Veilleux & Osterbrock (1987). The horizontal solid 
hues demarcate Sy2's from LINERS (Heckman 1980). Dotted hues are upper hmits from 
Dopita et al (2000) for zero age instantaneous bursts, dot-dashed hues are upper hmits for 
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Fig. 10.— Zoom of the STIS MAMA spectrum of F05 189-2524, showing the extended, 
asymmetric Lya hne profile (left) with the more compact, symmetric NV A1240 line. The 
vertical bar denotes 150pc. 
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42.0 I \ \ I \ \ \ \ \ \ \ \ \ \ \ \ \ I \ \ \ \ \ \ I \ \ \ \ \ \ \ \ \ \ \ \ 

120 130 140 150 

Wavelength (nm) 

Fig. 11. — Comparison of the UV spectrum of F01003-2238 with a Starburst99 instantaneous 
burst model with age 3Myr and an IMF slope of 2.0. The observed spectrum has been dered- 
dened using Equation 3 and is plotted in the rest-frame. The model has been offset below 
the spectrum for clarity. The gap in wavelength coverage for the model between llSOAand 
1200Aarises because Starburst99 does not produce a model spectrum in this wavelength 
range. 



-so- 




ns 



120 



125 130 135 
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Fig. 12. — Comparison of the UV spectrum of F15206+3342 knot 1 (top) with representative 
Starhurst99 spectral synthesis models. A 6Myr burst with an IMF slope of 2.35 (dashed line) 
is a good fit to the NV emission and a reasonable fit to the Lya absorption profile, but cannot 
explain the other absorption features. A 20Myr burst with an IMF slope of 3.3 (dot-dashed 
line) is a good fit to the Lya absorption, but cannot explain the NV emission and is a very 
poor match to the absorption features. Even combining the two models (solid line) does 
not produce a good fit. Note we are not invoking multiple star bursts in this spectrum, the 
models are shown simply for comparison. 



